Adiabatically steered open quantum systems: Master equation and optimal phase 
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We introduce an alternative way to derive the generalized form of the master equation recently 
presented in Ref. [J. P. Pekola et al, Phys. Rev. Lett. 105, 030401 (2010)] for an adiabatically 
steered two-level quantum system interacting with a Markovian environment. The original derivation 
employed the effective Hamiltonian in the adiabatic basis with the standard interaction picture 
approach but without the usual secular approximation. Our approach is based on utilizing a master 
equation for a non-steered system in the first super-adiabatic basis. It is potentially efficient in 
obtaining higher-order equations. Furthermore, we show how to select the phases of the adiabatic 
eigenstates to minimize the local adiabatic parameter and how this selection leads to states which 
are invariant under a local gauge change. We discuss also the effects of the adiabatic noncyclic 
geometric phase on the master equation. 
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I. INTRODUCTION 

The adiabatic theorem P, Q has been one of the 
workhorses of quantum physics for decades. It states that 
if the external control parameters of the system Hamilto- 
nian vary slowly enough, the system remains very accu- 
rately in one of its initial instantaneous eigenspaces. As 
slowly varying quantum systems appear in many fields of 
physics, a multitude of applications for the theorem ex- 
ists. In recent years, adiabatically steered quantum sys- 
tems have attracted a lot of interest due to their connec- 
tion to geometric phases in cyclic evolution [3H5J- These 
phases provide a potential alternative to quantum infor- 
mation processing in which the quantum gates are 
implemented by purely geometric means [9Ml8|. This 
geometric quantum computation has been shown to of- 
fer inherent robustness against control errors (lfjl - l^ due 
to the fact that geometric phases depend only on some 
global geometric properties. Different ways to describe 
geo metric phases in open systems have been introduced 
[23M28I I and methods to account for the effect of the en- 
vironment on the system evolution have been studied 
flfjl. I29l436| along with techniques to reduce the unwanted 
noise. However, a consistent description of the combined 
effect of adiabatic steering and noise was missing un- 
til recently, when a master equation was introduced in 
Refs. H3,s||. 

In the approach of Ref. [33|, it was shown that the 
typically applied secular approximation [32l . |39| is not 
suitable in describing adiabatic evolution. Taking into 
account all the relevant contributions leads to a mas- 
ter equation which ensures relaxation to a proper basis 
and shows that the ground-state dynamics are not influ- 
enced by zero-temperature Markovian noise in the adi- 
abatic limit. Thus the system exhibits inherent robust- 
ness. The master equation derived in Ref. [37j was gener- 
alized to hold for a g eneric system-environment coupling 
operator in Ref. [38fl . Furthermore, the master eq uation 
was applied to describe Cooper pair pumping [4fj-|42| in 



Refs. (37|,|3g. 

In this paper, we introduce an alternative derivation 
of the master equation for adiabatically steered quan- 
tum systems coupled to a Markovian environment. Our 
derivation is based on utilizing a non-steered master 
equation in the first super-adiabatic basis. We show that 
the master equation we obtain is the same as in Ref. (38| . 
Our method is potentially more efficient in obtaining 
higher order expansions in the adiabatic parameter. In 
addition, we introduce a way to select the complex phases 
of the adiabatic basis states such that the local adiabatic 
parameter is minimized leading to vanishing diagonal el- 
ements for the operator describing the steering. We show 
that this selection results in locally phase invariant basis 
states. Finally, we discuss how to account for the time- 
local accumulation speed of the geometric phase in the 
environment-induced transitions. 

The structure of this paper is as follows. In the next 
section, we introduce our model describing the open 
quantum system. In Sec. IIII1 we derive the master equa- 
tion for a non-steered system subject to decoherence. In 
Sec. IIV1 we use the non-steered master equation to ob- 
tain the full master equation for adiabatic steering. In 
Sec. [V] we introduce the optimal phase selection for the 
adiabatic eigenstates and demonstrate the main impli- 
cations of such a selection. We conclude the paper in 
Sec. EE 



II. MODEL 

We consider a quantum system with a Hamiltonian Hs 
which depends on a set of real control parameters {qk} 
that vary in time. The system is assumed to be interact- 
ing with the environment so that the total Hamiltonian 

is 



H(t)=H s (t) + V(t)+H E , 



(1) 



2 



where V(t) is the coupling between the system and its 
environment and He is the Hamiltonian of the environ- 
ment. We assume that the coupling is of the generic form 
V = A<E>X(t), where A is the system part of the coupling 
operator and X(t) acts in the Hilbert space of the envi- 
ronment. Let |m; q(t)) be the instantaneous eigenstate of 
Hs{t) and E m (t) the corresponding eigenenergy defined 
by Hs[q(t)] \m;q(t)} = E m [q(t)] \m;q(t)}. In the context 
of adiabatic evolution, {\m;q{t)}} is referred to as the 
adiabatic basis. We assume that the adiabatic states are 
normalized and non-degenerate. 

Let the Hamiltonian Hs(t) be diagonalized in a fixed 

basis {|t71/)} using the eigendecomposition as Hs(t) = 

&{t)H s (t)D(t), implying that (n f \H s (t)\m f ) = 
E m (t)5 nm . We define a similar transformation for the 
total density operator p(t) in the Schrodinger picture as 
p(t) = b^(t)p{t)t)(t). It follows from the Schrodinger 
equation that the evolution of p(t) is governed by the 
effective Hamiltonian for the adiabatic basis 



H (1 \t) = H s (t) + hw(t) + V(t) + H E 



(2) 



where V(t) = D^{t)V{t)D(t) = D**{t)AD(t) ® X(t) and 

w(t) = —iD' (t)D(t). Omitting the environment and as- 
suming adiabatic evolution, a more accurate approxima- 
tion for the exact evolving state is achieved if the adia- 
batic states are corrected by 



\Sm;q(t)) = -iH ^ \k;q{t)) - 



k^rn 



Ek 



(3) 



in the first order in the perturbation theory. The basis 
formed by the corrected states {\m) + \Sm)} is usually 
referred to as the first superadiabatic basis [3j| . 

We introduce the local adiabatic parameter as a(t) = 
H\\w(t)\\/A(t), where we compare the Hilbert-Schmidt 
norm of the operator arising from the adiabatic evolution 
||u)(f)|| = y/Ti s{w(t)^w(i)} to an instantaneous mini- 
mum energy gap in the spectrum A(i). Here Trs denotes 
the trace over the system degrees of freedom and in the 
following we will use Tr e to denote the trace over the en- 
vironment degrees of freedom. The parameter a(t) yields 
typically a good estimate for the degree of adiabaticity of 
the evolution (25l . l38| . In cyclic evolution with the period 
T, the parameter scales as 1/T and, thus, in adiabatic 
evolution we should have a(t) -C I • 



dard method [37H39j to obtain the master equation for a 
non-steered system as 



d Pg a 
dt 



= -(Tge + T eg )p gg + sRe{f p ge } + r e9 , (4) 



and 



- iUJQlPge - (f + + f -)p gg 



dt 



(5) 



where p rs = (r\ps\s) with r, s £ {g,e}, and u>o\ = (E e 
E g )/h. The transition rates are defined as 



r 



(e\A\g)\ 2 

(z\A\g)\ 2 



Sx(—uoi), 



r n = 



9 - h? 

e\A\g) ({g\A\g} - (e\A\e)) 



Sx(0), 



(g\A\e) ((e\A\e) - (g\A\g)) c 

1 ± = T^to OX(±^0l), 



(e\A\e) 



2h 2 

{g\A\g)\* (g\A\g) (e\A\e) 



2h 2 



r 



2ft 2 

(9\A\e) 
2ft 2 

(9\A\ef 
2ft 2 



ft 2 



•Sx(uoi), 



Sx(-^oi). 



The spectral density is denoted by Sx(w) = 
dTTr E {pEX(T)X(0)}e lulT . Note that we neglect the 
drive, i.e., omit all terms proportional to w. Further- 
more, Eqs. and ([5]) include all the nonsecular terms 
neglected in the usual application of the approach of 
Ref. [H. 

For details concerning the derivation of Eqs. (0| and 
([5]) see Appendix. Especially, we neglect the possible 
imaginary parts of the transition rates, i.e., the Lamb 
shift, by assuming that the system time scales are longer 
than the system autocorrelation time. 



IV. MASTER EQUATION FOR ADIABATIC 
STEERING 



III. MASTER EQUATION FOR A 
NON-STEERED SYSTEM 

Let us study the dynamics of a generic non-steered 
two-level quantum system coupled to its environment. 
Denote the ground and excited states of Hs in the 
Schrodinger picture as \g) and |e), respectively, with cor- 
responding eigenenergies E g and E e . We apply the stan- 



We aim to derive the the full master equation for the 
system coupled to its environment in adiabatic steering 
using the master equation for a non-steered system. De- 
fine a unitary transformation D\{t) making Hs(t) + hw(t) 
diagonal in the fixed basis {|0) , |1)}. Thus the evolution 
of the density matrix p^ = b\pbi = D\£>^ pbb\ is 
governed by the effective Hamiltonian for the first super- 



3 



adiabatic basis 



H {2 \t) = H {2) {t) + Mh(t) + V w {t) + H 



(2)/ 



(6) 



where H 
it 



(2) 



s (t) = D\(t)[H s (t) + hwWD^t), VW(t) = 

Z){(t)F(t)Di(t), and = -it>\ (t)j6 x (t). 

Assume that the super-adiabatic correction, u)i , is neg- 
ligible with respect to the adiabatic one so that we can 

write Eq. (JTSJ) as ^(f) « ^ 2) (i) + ^ 2 )(f)+ff s . Since 
this Hamiltonian describes effectively a non-steered sys- 
tem, we can employ the approach of Sec. IIIII to write a 
master equation similar to Eqs. (fj| and ([5]) as 



j (2) 

dt 



r (2)) (2) 



(2) (2)| +r (2) 
Fpe J * eg ) 



(7) 



sume that the system is driven adiabatically but does 
not necessarily remain in the ground state at all times, 
we cannot assume the density matrix elements p ge to be- 
come small enough to be neglected due to their order in 
this limit. Hence, we are only considering a(t) as a small 
parameter and neglect all terms with a 2 or higher order. 



If we rewrite p gg and p ge using Eqs. ([9]) and (|10|) . the 
master equation becomes 



Pgg - 2" 



woi 

f r (2h 

~ eg ) 



,(2) 



Pgg 



^(w* p ge ) 



.w 



g> 



2-f- Pgg 



Woi 
OJQl 



r( 2 ) 

eg 



(11) 



and 



a (2) 



• (2) , 
^01 Pge 

<2) 



- (f (2) 



fL 2) )4 2) 



,(2) 



, , p(2) 

2 2 v 

+ ' 



Pge 



(r 



(2) 



(8) 



where we have marked the relevant terms in the super- 
adiabatic basis with the superscript ^ to avoid confusing 
them with the adiabatic ones. The transformation D\(t) 
can be approximated using the perturbation theory for 
the adiabatic correction. This results in the superadia- 
batic eigenstates which we can obtain from Eq. (|3|) up to 
the linear order in a(t) in the two-state model as 



l.9 (2) ) = Iff) - \e) ^ 



and 



,(2) 



E„ 



- hvjgg and E, 



(2) 



(9) 



(10) 



E P 



with the eigenenergies Eg — 
hw ee , respectively. Here, we denote the matrix elements 
of the adiabatic correction as w sr = — i (s\r), where r, s £ 
{g,e}. Thus, the super-adiabatic energy gap up to this 

(2) 

order becomes Wqi = Woi + ( w ee — w gg ). 

The matrix elements in Eqs. ([7]) and (J5| can be writ- 
ten using the super-adiabatic eigenstates to obtain the 
master equation for adiabatic steering up to the linear 
order in a(t). We restrict our derivation to the linear 
order since in the adiabatic limit, a(t) — > making the 
contributions beyond the linear one negligible. If we as- 



and 

Pge 



^Pgg 

i[u l + (Wee 



(ff 4 

n(2) 
L eg 

2 

(r<?> 4 



f(2)> 



2f^Pgg 



Woi 



Pgg 



woi 



.(2) 

ge 
2 



r (2) 



,(2) 



(^Pge + 


2^Pgg 


_ «v 




Woi 


w ge 

—Pgg 


_ 




WOl 


Woi / 



(12) 



where we have neglected all terms of order a 2 or higher, 
except in the terms which treat further below. We 
can solve p gg and p ge from these equations to obtain the 
full master equation. In addition, we employ Eqs. Q and 
(jlOP to rewrite the rates in the super-adiabatic approx- 
imation. To present the full master equation, we adopt 
a notation which will not reduce the generality of the 
equations but simplify them. In the nested commutator 
expression for the master equation for a non-steered sys- 
tem [see Eq. (|A.1[) in Appendix], the coupling operator 
is only found in places where it is commuting with other 
operators and, hence, provided that the Lamb shift is ne- 
glected, we can add any operator comparable to the iden- 
tity operator to it without affecting the nested expres- 
sion. Thus, the system part of the coupling operator can 
be choosen traceless in the two-state basis. We will adopt 
this convention by introducing mi = (g|^4|<7) = — (e|A|e) 
and 777-2 = (<?|^4|e). Notice that mi € R whereas € C 
in the case of a general coupling operator. The master 
equation up to the linear order in a(t) and the quadratic 
order in the system-environment coupling becomes 



4 



p gg = -2^sm(w* e pg e ) + 5 ( uj a 1 ) | m 2 1 2 - [S(—uoi) + S(ujoi)]\m 2 \ 2 pgg + 2[5m(m 2 )Sm(p S e) + 5fte(m 2 )3fie(p ge )]5(0)mi 

- 2 2g(0) ~ giz^l) ~ g( " 0l) {[Sm(m 2 )^m( % ) + SefmjJSe^)] [3m(m 2 )3m(/), e ) + 5Re(m 2 )Ke(p 9e )]} 

£*>oi 

25(0) — 5(— cjoi) — S(woi) , . . . \ I cv\ / \ <x\ / 
+ 2 { L ~jm(m2)<sm(Wg e ) + !Re(m 2 )Ke(w ffe )}mip gg 

— 2 - ^ ^ mi{^sm(m2)^'m(wg e ) + 3Re(m2)5Re(w se )} 

(13) 



and 



f)g e = iWg e {2pgg ~ l) + « (l^e ~ Wgg)p ge + iOJ 01 Pg e ~ 5 ( W() 1 ) TTl 1 TO 2 + [5(-W l) + 5 ( ^0 1 ) ] ^ 1 »™2 P BB ~ 2S(0)m\p ge 

- i[S{-u) i) + SVoi)]«i2[3m(p 9e )Ke(m 2 ) - Ss(m 2 )$te(p ge )] 
O 2g(0)-g(-o; 01 )-g(a;o 1 ) 2 g(0) - SKQ 2 

- 2 myWgePgg + 2 myWge 

CJoi Woi 

5(— Woi) — 5(cJoi) rc v. / ncts / s \<tw M 

25(0) - 5(-woi) - 5(woi) r . rc x. / w , , , > 0/ \ i 
m 1 {irn2[-*sm(Wg e yKe(p ge ) - \sm{pg e yKe(w ge )\ 



woi 

[3m(m 2 )3m(w se ) + Ke(m 2 )5ie(w; ge )]p ge }. 



(14) 



Here, we applied a shortened notation for the spectral 
densities 5(w) = Sx(w)/h 2 . We would like to empha- 
size that in this section, we assume that the system is 
externally steered, i.e., the system Hamiltonian is time- 
dependent. Even though we do not explicitly make the 
approximation of adiabatic rates (39| . i.e., assume the 
evolution time is much longer than the environment au- 
tocorrelation time so that woi, mi, m 2 , and the matrix 
elements of w vary slowly in time, the approximation is 
implicitly assumed. This assumption stems from the fact 
that we use the master equation for a non-steered system 
in the linear order in a{t). Furthermore, we have ne- 
glected corrections proportional to d u S(uj)\u=±u> 01 when 
estimating the power spectra above since they correspond 
to driving induced Lamb shift contributions which are ne- 
glected also in the derivation in Refs. (37l.[38|. In Sec. fVl 
we show that w gg and w ee vanish from Eq. (|14[) for a 
specific choice of the phases of the adiabatic basis states. 

Let us explicitly reformulate and assess the range of 
validity of the master equation [see Eqs. P^|) and (fl"4"|) ]. 
The approximation of adiabatic rates requires that the 
external drive does not change the Hamiltonian govern- 
ing the system on the time scale of the memory time of 



the bath, i.e., we assume that r corr <C l/||w)|| where r corr 
is the environment correlation time. Furthermore, we 
require that the dynamics of the density opertor occur 
on time scales longer than the environment autocorrela- 
tion time so that t co1t <C 1/(jl>o, where l/u>o is a typical 
system transition time relating to the off-diagonal ele- 
ments of the density matrix. As mentioned in Appendix, 
we apply the Markov approximation by assuming that 
Tcorr < I/7 where I/7 is a typical relaxation time of the 
system. Moreover, the assumptions we employ demand 
that l/u>o l/7i l/HHI- Combining the used approxi- 
mations leads to a requirement of time scale separation 

<c i/ujq < 1/7, i/|H|. 

Remarkably, our master equation is identical to that 
derived in Ref. j38| . however, the manner in which the 
master equation was derived is different. In Ref. (38| . 
one starts from the effective Hamiltonian for the adia- 
batic basis presented in Eq. ([2]) and formulates a nested 
commutator expression for the derivative of the reduced 
system density operator in the adiabatic basis applying 

hw(t) + V(t) as the perturbation 



5 



d ^- = i[h(t),w I (t)]-^ E \l Q 




h{t)AMt')Mt")]] Mt) 



(15) 



r 



in the interaction picture. Using this operator directly, 
results in the same master equation to the one we ob- 
tained. Thus, we find that with respect to adiabatic 
temporal evolution it makes no difference whether one 
uses the effective Hamiltonian for the adiabatic basis and 
takes hw(t) + V(t) as the perturbation as was done in 
Refs. [13, or whether one uses our approach to express 
the effective Hamiltonian for the super-adiabatic basis as- 
suming that the super-adiabatic correction is small, thus 

taking V^ft) as the perturbation and writing the super- 
adiabatic basis states up to the linear order in a(t). Our 
discovery reaffirms that the super-adiabatic basis approx- 
imates the exact evolving state in the next order in a so 
that using only the bath coupling as the perturbation will 
result in describing the dynamics in the same order as the 
effective Hamiltonian for the adiabatic basis does. This 
result is an important consistency check for the master 
equation derived in Refs. [13, [H|, see Eqs. $F5§ and ([Til). 

The original way [13, H3 of deriving the full master 
equation can be extended to obtain master equations in 
higher orders in a by applying the nesting procedure it- 
eratively [see Eq. ([T5j) [ . Our method can be used as well 
to obtain higher-order master equations by using higher- 
order perturbation theory to write the matrix elements 
required to utilize the master equation for non-steered 
systems. Since our technique is based on applying alge- 
braic operations, it is potentially simpler to obtain higher 
order equations with it than with the nesting procedure 
which results in complicated integral expressions. 

As an example, let us illustrate how to obtain the mas- 
ter equation up to the second order in a(t) using our 
method for a two-level system. We begin by defining 

~ (2) 

a unitary transformation D 2 it) making H s (t) + hw\ (t) 
diagonal in the fixed basis {|0} , |1)}. Thus the relevant 
density matrix becomes = D\p^D 2 and its evo- 
lution is governed by the effective Hamiltonian for the 
second super-adiabatic basis 



#( 3 )(t) = Hf\t) + hw 2 (t) + V^{t) + H E , 



(16) 



•e H ( s\t) 



b\{t){Hf\t) 



Dl(t)V^(t)D 2 (t), and w 2 = ~ib\{t)b 2 (t). Assuming 
that w 2 is negligible with respect to wi , we can again re- 
sort to the approach of Sec. IHII and write a master equa- 
tion similar to Eqs. ([4]) and ([5]). Finding the relevant ma- 
trix elements is slightly more complicated than in the case 
of the master equation up to the linear order in a(t). We 



essentially wish to describe D(t)Di(t)D 2 (t) \m), where 
\m) is a fixed state, up to the second order in a(t). We 
first write the first super-adiabatic eigenstates up to the 
second order in a(t) as 



and 



9 {3) ) = \ 9 ) - \e) ^ 



| e <B)) = | e> + \ g) ^1 + {g) 



WeeW ge 



woi 



01 



WggWge 



-01 



(17) 



(18) 



with the eigenenergies E g 3 ^ = E g + hwgg — fr\w ge \ 2 /luqi 

and — E e + hw ee + h\w ge \ 2 /woi> respectively. Ad- 
ditionally, we must account for the lowest order in 
||t/>i(£)|| ~ a(t) 2 to obtain the states of interest up to 
the second order in a(t) as 



and 



£>(i)A(i)£> 2 (i)|l>«|e( 3 )} + | ff ) 



W01 



WOl 



(19) 



(20) 



and the energy gap as uj^ = l 
{w x> ee ~ wi.gg). Here we 
it 



(2) 
01 



2k ge | 2 



7^01 + 
denote Wi. gg = 

-i(0|£>I(i)Z)i(i)|0>, w hee = -i(l|^(t)A(t)|l> and 

w i,ge = —i (0\D\(t)Di(t)\l) . Using these definitions, one 
could write the matrix elements in the relevant master 
equation for a non-steered system to obtain the master 
equation for adiabatic steering up to the second order in 
a(t). 



V. OPTIMAL PHASE SELECTION 

Assume that we have nondegenerate adiabatic eigen- 
states \g) and |e) that are normalized and smooth during 
the temporal evolution, and that they can be obtained 
from the fixed states with a unitary transformation as 
D |0) = \g) and D\l) = \e). Thus, the operator deter- 
mining the adiabatic evolution becomes w = —iD^D [see 
Eq. J2J], However, the choice for the complex phases of 
the states is essentially arbitrary, and hence its effect on 
w should be studied. This means that we could also work 



(i 



with basis states that differ from \g) and |e) by phase fac- 
tors which depend on the point in the control cycle. In 
this section, we show that this degree of freedom can be 
used to minimize the local adiabatic parameter, while in 
the corresponding master equation it leads only to renor- 
malized matrix elements. Furthermore, we discuss how 
to account for the noncyclic geometric phase leading to 
a slightly different but generally more feasible phase se- 
lection. 



Let us choose new phases for the states by multiplying 
them with phase factors e tXg and e lX ° , where X g , X e 6 M, 
so that a phase selection operator is defined as 



in Eq. 



to vanish, yielding 



h = e lXg |0) (0\+e iX ° 



(21) 



yielding the new transformation as D = DVl. Notice 
that the new states defined by the transformation are 
still eigenstates of the original Hamiltonian. With this 
transformation, the operator for the drive becomes 

w = -ib^D 

= -iCtf b^ \bn + bh) ( 22 ) 
= \ g |o) (o| + x e |i> (i| - ih^b^M, 

where we have used the unitarity of D. The matrix ele- 
ments in the phase shifted basis become 

it,' 



<0|W|0) = X g - l (0|£>T£>|0) = Xg + Wgg, 

(l\w\l) =X e -l (\\tflb\\) =X e + Wee, 

(0\w\l) = -ie l ^- x ^ (OlDt^ll) = e l ^- x ^wge, 
(l\w\0) = -ie^ x o- x ^ (l\D^b\0) = e l ^- x ^weg. 



(23) 



Thus, the phase shift induces a shift in the diagonal ele- 
ments and a phase shift in the off-diagonal elements. 



Since the choice of the phases of the eigenstates 
changes the w terms of the master equation and affects 
the quantum evolution, it has to be fixed by physical 
reasoning. One way to avoid artefact effects arising from 
the phase choice is to minimize the Hilbert-Schmidt norm 

|w|| = <J Trgjiytii)} at each time instant. For this task, 
it suffices to minimize 



Tc s {&w} = | (OHO) | 2 

+ l(iHo>| 2 - 



(iMi)l 2 
:o|^|i)l 2 . 



(24) 



The last two terms consist of the off-diagonal terms and 
the phase selection has no effect on them. Thus, the 
minimum is found when we select the diagonal elements 



\{t) 



I" Jt'w gg (t') + X° g =i [ dt' (g\g)+X gl 

Z Z (25) 

Xe(t) = - / dt'w ee (t') + X° e = i dt' (e|e) + A°, 
Jo Jo 

and 



h \ = e i(K-K) e i So dt'[w gg (t')-w ee (t')] t 



'ge 



e ^°-K) e ^ dt '^\g)-(e\e)] w _ ^ 



The absolute phases are not fixed since we have a de- 
gree of freedom in the selection of the constant parts 
A° and X ( g . Notice that the primary selection of the 
smooth eigenstates \g) and |e) determines the accumulat- 
ing phase. We denote the integrals as simply over time, 
but one should bear in mind that they contribute a path 
in the parameter space. Used in conjunction with our 
master equation, the above phase selection results in w gg 
and Wee vanishing in Eq. (|14p . Furthermore, it minimizes 
the local adiabatic parameter a(t) = ||tD(t)||/o;oi(£). 
Thus we refer to it as optimal phase selection although it 
may not yield the most accurate evolution as discussed 
below. 

Utilizing the optimal phase selection scheme with our 
master equation requires a careful consideration of the 
used approximations, in particular, the approximation 
of adiabatic rates. We used this approximation in the 
derivation of the master equation [see Eqs. (fT3")) and (fl"4"])] 
to state that cjoi> mi, 7712, and the matrix elements 
of w vary slowly in time. With the optimal phase se- 
lection, the corresponding parameters are (Dqi = ujqi- 



mi 



a i(X e -X g ) 



Wgg 



0, and 



mi, m 2 

Wge = e l ( x '~ Xa ^Wg e . Since the approximation of adia- 
batic rates applies for the derivatives of the accumulating 
phases defined in Eq. (|25|). any possible shift in the power 
spectra induced by the optimal phase selection can be ne- 
glected since it would only lead to higher order terms in 
the master equation. Thus, the master equation can be 
directly used by replacing the original variables with the 
phase shifted ones. 

With the optimal selection, the phase shifted ba- 
sis states become \g) — e lX ne~ fo dt (sis) \g) and |e) = 

e i\° e -f*dt'{e\e) Thus we haye _ ^g|g^ _ q in _ 

dependent of \g) and |e) and the selection renders the 
phase shifted states to be also invariant under a local 
gauge change, i.e., \g) — > e 1 ^^ 1 \g) has no effect on \g) 
and |e) —¥ e" 1 ^ |e) has no effect on |e) where j5{t) and 
r](t) are any smooth functions. For a closed path in the 
parameter space 7, we have 



K(tb) ~ X g (t a ) = i <b (g\g), 

J y 

X e {tb) - X e {t a ) = 1 f (e\e) , 

J 7 



(27) 



7 



where we have denoted t a and % as the virtual start- 
ing and ending time instants for the path, respectively. 
These are the Berry phases accumulated over the path 
for the phase shifted basis states and, as such, cannot 
be removed by any continuous local gauge change Q. 
Thus, selecting the optimal local phase for a closed loop 
in the parameter space implies a gauge- invariant accumu- 
lated phase at the end of the loop. However, the optimal 
phase selection neglects the effects of the accumulation 
speed of the geometric phase on the environment-induced 
transitions in the master equation. Although this effect 
is negligible in the adiabatic limit, it is useful to provide 
means to take it into account consistently. 

To go beyond the optimal phase selection, we need to 
study the adiabatic noncyclic geometric phase acquired 
during the evolution f43l . [44j. The adiabatic noncyclic 
geometric phase acquired by the nth eigenstate [44| can 
be defined as 

7„(t) = axg{{n;q\t )\n;q{t)}} 

r* ( 28 ) 

+ 1 / dr (n; q(T)\d T \n; q(r)) , 
Jo 

where |n; q{to)} is an essentially arbitrary reference eigen- 
vector not orthogonal to \n;q(t)}. The phase 7„(i) is 
invariant under any local phase transformation of the 
basis states and independent of the speed, at which we 
traverse the open path in the control parameter space. 
To allow for the accumulating phase of the basis vec- 
tors to describe the noncyclic geometric phase, one can 
define a phase transformation of the eigenstates |n) = 

e -i a rg{<n';g(to)|n';«(t)>} where n g e } and ^ de _ 

scribes an eigenstate with an arbitrary continuous local 
phase. This transformation produces 

w nn = -d t a,ig{(n';q(t )\n';q(t))} -i(n'\n'} 
= -7n(*)- 

Thus the time-local accumulation speed of the geomet- 
ric phase can be made to appear in the master equation 
[Eqs. (fT5|) and (TH)) ]. However, the accumulation speed 
in Eq. (f2T))) is dependent on the reference point \n'; q(to)) 
and, thus, redefining the reference changes typically the 
accumulation speed at all times. This can become a prob- 
lem since the accumulation speed affects in general the 
resulting physical quantities which should not depend on 
the choice of the reference point. A possible way to cor- 
rect for this inconsistency is to use the geodesic approach 



VI. CONCLUSIONS 

We deviced a way to derive the full master equation 
for adiabatically steered quantum systems in the two- 
state approximation under the influence of decoherence 
starting from an interaction-picture-based derivation, in 
which the external drive was first omitted. The full mas- 



ter equation was obtained by approximating the trans- 
formation to the superadiabatic basis using the pertur- 
bation theory and exploiting the master equation for the 
non-steered system. We showed that the master equa- 
tion we obtain this way is the same as the one obtained 
in Ref. [38| by a longer calculation. We concluded that 
our manner of obtaining the master equation is a con- 
sequence of the superadiabatic basis approximating the 
exact evolving state in the linear order in the adiabatic 
parameter a(t). Furthermore, there is no need to eval- 
uate high-order nested commutators of integrals in our 
method if it is extended beyond the linear order in a(t) 
as opposed to the method in Refs. [13, HI]- A detailed 
study of the efficiency of these two approaches is left for 
future research. 

There exists a gauge degree of freedom in the choice 
of the phases of the basis states during the evolution. 
We demonstrated a way to choose the phases in a way 
which minimizes the local adiabatic parameter and sim- 
plifies the derived master equation. We showed that this 
choice produces basis states which are invariant under a 
local gauge change. Furthermore, we discussed how to 
account in a gauge invariant manner for the effects of the 
accumulation speed of the adiabatic noncyclic geometric 
phase on the environment-induced transitions. 
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Appendix: Non-steered master equation in the 
two-state basis 

The reduced system density matrix in the interac- 
tion picture <?i(t) = Trs{/5/(t)} can be used to de- 
rive the relevant master equation assuming a station- 
ary environment, i.e., -ffi = j:[pe,He] = 0. We de- 
fine the operators in the interaction picture as Zi{t) = 
e i6Bt / h Ul(t,0)Z(t)U s (t,0)e- iAl!!t / H , where Z(t) is the 
operator in the Schrodinger picture and Us(t,0) is the 
time-evolution operator. For a time dependent system 
Hamiltonian, the time-evolution operator is Us(t,0) = 
7- e -i/o^s(T)dr/ft but simp i ines t0 Jj s (t,0) = e - ikst ' h 
for non-steered systems studied in this Appendix. 

If we assume that the system interacts weakly with the 
environment, the master equation acquires the standard 



s 



form (Redfield equation j45|) 



d&i(t) 
dt 



h 2 



dt'Ti- E ([[aj(t) ® p E ,Vi(t% Vr(*)]). 

(A.l) 

in the interaction picture, where we have utilized the 
Born-Markov approximation [45j . The transformation 
from the interaction picture to the Schrodinger picture 
unfolds when we employ 



Derivating Eq. (|A.2|) yields the transformation of the 
derivative as 

^ = ife(*),^(*)]+^(t,o)^^(t l o). 

(A.4) 

Using Eqs. (|A.1[) and (|A.4[) . we define the diagonal matrix 
element 



Ps(t) = U s (t,0)a I (t)Ul(t,0) 



(A.2) 



which can be used to obtain the density matrix transfor- 
mation componentwise as 



Pgg(t) = Vl,gg(t), 
Pee(t) = 0- I>ee (t), 

Pgeit) = e iumt a I>ge (t), 



(A.3) 



t dp s (t) 
<J>\-&-\9) = 



h 2 



dt'Tr E {(g\[[cri (t) ® fa, Vj (*')] , V> (t)] }, 



(A.5) 



and the off-diagonal matrix element 



iuJmp ge (t) 



iojQit 



¥ I dt'Ti E {{g\[[a I {t)®p B Mt%Vi{t)]\e)}, 



(A.6) 



r 



for the non-steered master equation. Notice that our 
derivation is based on assuming that the system relax- 
ation time is long compared to the environment correla- 
tion time T cc ,rr so that the environment has no memory, 
i.e., we are in the Markov regime [3{|. This allows us to 
neglect any variation of 07 (t) between times t and t+T corr . 
The integral expressions in Eqs. (|A.5|) and (|A.6[) simplify 
to give Eqs. and ([5]) when we expand the commuta- 



tors, use the closure relation for the adiabatic basis, and 
utilize Tv E {p E X(t')X(t)} = ^S x (uj)e-^ (*'"*) and 
Eq. (|A.3[) . Furthermore, we assume that the system time 
scales are longer than the system autocorrelation time to 
approximate the spectral densities in the remaining inte- 
gral expressions. This treatment leads to neglecting the 
Lamb shift. 
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